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Using the Gell-Mann-Oakes-Renner (GOR) relation and 
semi-classical arguments, we show that the bulk quark spec- 
trum in QCD exhibits a variety of regimes including the er- 
godic one described by random matrix theory. We analyze 
the quark spectral form-factor in the diffusive and ballistic 
regime. We suggest that a class of chiral transitions in QCD 
is possibly of the metal-insulator type, with a universal spec- 
tral statistics at the mobility edge. 
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1. The chiral properties of QCD in a small Euclidean 
volume are constrained by the way chiral symmetry is 
spontaneously broken in the infinite volume limit [jj]. 
This observation carries to the quark spectrum averaged 
over gauge configurations of fixed winding number, in the 
form of microscopic sum rules . The sum rules can be 
rederived using chiral random matrix models ||. The 
microscopic sum rules reflect on the universal properties 
of the eigenvalue distribution and correlations a quantum 
leap apart Q. 

The interesting results achieved by random matrix 
models in describing the microscopic regime of chiral 
QCD, prompt us to ask whether we could think about 
other generic or universal aspects of the quark spectrum, 
aside from that of random matrix theory. In this let- 
ter, we would like to provide some new insights to the 
microscopic regime using semi-classical arguments. We 
will use these new insights to analyze the quark spectral 
form-factor in the diffusive and ballistic regimes (non- 
universal). We then argue that a number of chiral tran- 
sitions in QCD are possibly of the metal-insulator type, 
and suggest that at the edge of the metal-insulator tran- 
sition a universal spectral statistics is expected. 

2. To be able to organize the various energy scales of 
interacting quarks in a finite Euclidean volume V, we call 
upon the two-flavor GOR relation H, 



F 2 = (qq) = — irp(0) , 

mi mt 



(1) 



where F is the pion decay constant, m T the pion mass 
and (qq) = — S is the quark condensate. The last identity 
in (Q) follows from the Banks-Casher relation with 
p(0) the quark density of eigenvalues per unit four volume 
V = L 4 at zero virtuality (A ss 0). Recall that in a metal, 
the Kubo formula for the dc-conductivity a in terms of 



the diffusion constant D and the density of states pf 
at the Fermi level is (in units of e 2 /%) : a — Dpp . By 
analogy, we see that in QCD, D w _F 2 /E plays the role of 
the diffusion constant while the Fermi level corresponds 
to the zero virtuality point. This physical analogy will 
be confirmed below by a direct calculation. 

The appearance of a diffusion constant D allows us 
to organize the various stages of the disordered phase. 
For a quark with proper time t ~ 1/A, where A is its 
virtuality, the relevant time scales are: the Heisenberg 
time, tu = 1/A with A = 1/pV the typical quantum 
spacing; the ergodic time, T crg = 1/E C = L 2 /D, with 
E c the Thoulcss energy [Q; and the diffusion time = 
1/2toq with rag the constituent quark mass. For fixed 
V, the Ohmic conductance is a L = E c /A « F 2 L 2 . The 
relation D = VpT^/i suggests that the virtual velocity of 
diffusive quarks at A « is about vf ~ F(m,Q /S) 1 ^ 2 ~ 1. 
Diffusive quarks are the analogue of constituent quarks. 

We will refer to the quarks in the regime t > tjj as 
quantum, those in the regime r crg < t < ta as ergodic, 
those in the regime < t < r Glg as diffusive, and finally 
those in the regime t < Td as ballistic. We note that 
the time scales are ordered as tn ^> T crg ^> or V ^> 
y/V 3> 1 in units where F = E = mq = 1. The virtual 
quark spectrum in QCD bears much resemblance to the 
real electronic spectrum in disordered metallic grains [g| . 

3. An important tool for studying the quark level cor- 
relations in the QCD vacuum is the averaged two-level 
correlation function [01 



R A , E (s) = R{s) = 4r(KA+)KA-)) 



(2) 



measuring correlation between two eigenvalues A± in the 
virtuality window around A with width E, being apart 
from each other by As. v\ — 1/AF is the mean density 
in the band E. The averaging in (^) is over the gauge con- 
figurations A with the QCD measure. The quark eigen- 
values in a fixed gauge configuration A are denoted by 
%f[A]q n — \ n [A]q n , and the unaveraged density of states 
is 
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X k [A}). 



(3) 



In particular p(0) = (i/(0)). Eq. (g) is related to the two- 
point connected correlation function discussed in fll||i"c|] 
for constant It satisfies the sum rule J dsR(s)ds = 
and its Fourier transform is the spectral form-factor 
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Following the general semi-classical arguments by 
Gutzwiller [pH , Berry |l^] and others |l3|j , we suggest 
that for QCD 



K(t) 



2\t\A 2 
(2n) 2 f3 



p(t). 



(5) 



where p(t) is the quark return probability for fixed proper 
time t. In QCD with N c > 3 (number of colors) the 
quarks are in the (complex) fundamental representation. 
Their evolution operator along the proper time does not 
enjoy time-reversal invariance in the presence of chromo- 
magnetic interactions. In the semi-classical approxima- 
tion leading to (|j) this corresponds to f3 = 2 (only the 
diagonal terms in the closed orbit indices are retained). 
In QCD with N c = 2 the quark orbits are paired by the 
conjugation operator Ct-iK where C is charge conju- 
gation and K complex conjugation, hence (3=1 (diago- 
nal and off-diagonal terms within the pair retained). In 
QCD with Nf adjoint and real Majorana quarks, the or- 
bits are paired by CK || . If only one pair of orbits is re- 
tained (square root of the theory) then S = 7rp(0)/2 and 
/3 = 4. A similar observation applies to Kogut-Susskind 
quarks on the lattice. 

The normalized return probability p(t) of a quark from 
x(0) back to x(t) for fixed proper time t in a four volume 
V, is 



-*(Ai— A a )|*| 



V 2 fdXidX^ 

(Tr(S(x,y;zi)S\x,y;z 2 ))) (6) 
with Zi 2 — m — 1X1.2, and 



S(x,y;z) = (x\ 
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«y + iz 
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where N is the number of quark states in V. A proper 
regularization of the short paths is assumed. Setting 
A1.2 = A ± A/2 with A ps (near the zero virtuality 
point) then z\ = z — m — iA/2 and Z2 — z* . In this 
case, the quark propagator satisfies 

S*(x,y;z*) =-j 5 S(y,x;z)j 5 (8) 

and (Q) is seen to relate to the pion correlation function 

l Qb C«(x, y; z) = (Tr (S(x, y; z)i l5 r a S(y, x; z)i lb r b )) (9) 

For z — m, pion-pole dominance (long paths) yields 



(10) 



where the sum is over the pion momenta — n^ir/L 
in a symmetric box V — L 4 . Using the GOR relation (jl|) 
and the analytical continuation m —> m — iX/2, we find 



C n (x,y;z) fa — ^ 



JQ-{x-y) _ 



2|(gg)| 

-iX + 2m + DQ 2 



(11) 



with the diffusion constant D = 2F 2 /Y>. Inserting ( pd| ) 
into (||), and noting that J dA — E with E/A — N and 
p = l/AV, with E = \(qq)\ = irp, we conclude after a 
contour integration that 



p(t) =e -2-l*l^ e -^Q 2 KI 



(12) 



The prefactor e 2m|t| j g ex p ec t ec j damping for zero 
mode quarks. Heavy quarks with bare masses 2m > 
D(2ir/L) 2 (symmetric box) do not diffuse. 

The present arguments will allow us to analyze the 
spectral rigidity 



S 2 (7V) 



ds (A - \s\)R(s) 



(13) 



-N 



This is the variance in the number of virtual quark levels 
in an energy band E = NA, averaged over the QCD 
gauge configurations. Its dependence on N or E changes 
qualitatively with the virtuality band considered as we 
now discuss. 

4. In the ergodic regime r clg < t < tn the result ( p^ ) 
is dominated by the constant mode n — 0. Hence the 
return probability is p(t) — e~ 2m l'l, and 



R(s) 



1 



Ptt 2 (a 2 + s 2 ) 2 



(14) 



with a = 2m/ A. The corresponding spectral rigidity is 

'N 2 +a 2 ' 



= In 
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(15) 



for N ^> 1. For quark masses in the range a ~ 1 or < 
rnSy < 1 (V ~ iV) [E]J§| , the spectral rigidity is universal, 
S 2 (A) = 2lnN/(f3n]\ This is in agreement with the 
result of standard p5[ and chiral random matrix 

theory. For quark masses m ~ 1 in units where i 7, = S = 
1 , the spectral rigidity (|l^) deviates from random matrix 
theory. In QCD this mass range is already sensitive to 
strangeness. 

The result is also valid for quenched QCD as ( ]To[jTl| ) 
and the GOR relation hold in this case as well. For a ~ 1, 



the asymptotic of (14) is enough [[15 17| to show that the 
tail of the bulk level spacing distribution follows that of a 
Wigner surmise |15[ . The arguments in |L7| can be easily 
extended to general a. 

The present arguments only yield the perturbative part 
( |l4| ) of the two- level correlation function (||). The oscil- 
latory part well-known from random matrix theory [jl5| , 
requires a reassessment of (j|) including the short paths 
which we have not done in this work. In this context 
an interesting observation was made recently by Agam, 



2 



Altshuler and Andreev |Q, relating not only the per- 
turbative part of @ but also the oscillatory part to the 
spectrum of the diffusion operator (in our case (pi])), and 
more generally to the Perron-Frobenius operator. A simi- 
lar observation should hold in the context of QCD, where 
the short paths are just controlled by QCD perturbation 
theory, thereby establishing the chaotic character of the 
quark spectrum in QCD from first principles. 

Finally, we note that the relation (|^) in combination 
with (|12|) allows us to write 



1 



(sA + iDQ 2 + i2m) 2 



(16) 



in agreement with the perturbative result derived by 
Altshuler and Shklovskii jl9| for disordered electrons in 
metallic grains. In our case, the squared denominator in 
( filf ) follows from the exchange of 2-pions in the double 
ring diagram corresponding to the density-density corre- 
lation function (j^). This diagram is forced by G-parity 
and is dominant at large t. So QCD with N c > 3 cor- 
responds to the exchange of diffusons |l9[] , hence (3 = 2. 
QCD with N c = 2 and Nf > 2 allows for massless pions 
and baryons. In this case, the baryons are the analogue 
of the cooperons |]l9| , hence /3 = 1. 

6. In the diffusive regime Td < t < T org the result fll2] ) 
receives contribution from all modes in the box V. The 
outcome is still a classically diffusive motion of the quark 
in d=4, with a return probability 



P(t) 



-2m\t\ _ 



V 



{AllDtf 



(17) 



Setting 1/m > T crg that is T,/V < m < D/W and 
dialing the virtuality A such that 1/T crg < A < l/r<j will 
allow us to probe the diffusive part of the quark spectrum 
in QCD. In particular, the spectral form- factor (|5|) can 
be readily found. Its Fourier transform is seen to diverge 
for short times. Setting a short time cutoff at Td (the 
diffusion time), we obtain 



R(s) 



Z^(C + MAr d | S |)) 



{2n) 4 D 2 (3 



(18) 



where C = 0.577 is Euler's constant. We note that Ar^ = 
Td/tn 2> 1. In the diffusive regime the spectral rigidity 
is 



E 2 (A): 



VY 2 



4f3(2ir) 4 F 4 



(3 - C - ln(AAr d )) (NA) 2 (19) 



where N ~ E / A ^> lis the mean number of states in the 
energy band E. In the diffusive regime Erd < 1 so the 
spectral rigidity is always positive. However, it decreases 
with increasing N or E — NA. Although not universal, 
this result should be specific to the QCD quark spectrum 
in the quoted regime, and through the l/F 4 behavior 
amenable to chiral power counting M . 



7. In the ballistic regime t < Td the classical trajecto- 
ries of the quarks are shorter than the typical mean free 
path I = vpTd — vp/2mQ, which is about the constituent 
quark Compton wavelength. Although in this regime the 
density of states vary with the virtuality A, we will for 
simplicity ignore the variations in our case |2l| . For short 
times, the return probability is about constant but small, 
say p(t) ~ p + 0(t 2 /rj). Hence, 



R(s) 



/3(tts) 2 



1-0 
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(sATdf 



(20) 



A direct calculation of p may be averted by noticing that 
(O) implies 



dY 2 {N) 
dN 



N 



dsR(s) 



(21) 



—N 



for a fixed N = E/A, in the ballistic regime. Since (g) 
obeys the null sum rule on the entire support of eigenval- 
ues as required by the conservation of the number of en- 
ergy levels, it follows that the contribution from the bal- 
listic regime balances the contribution from the diffusive 
regime for L ^S> I, thereby fixing p. Setting the cutoff be- 
tween the ballistic and diffusive regimes at a; = D/ (I 2 A), 
we obtain 122JI 



R(s) = - 



V 



I 2 



1 



(22) 



We note that l 2 /(Drd) ~ 4, and the prefactor yields p « 
V/l 4 > 0, which is expected from free phase space in 
d = 4. In the virtuality range A ^> 1/t^ ~ rrtQ, the 
spectral rigidity for fixed E/ A = N ^> 1, is 



E 2 (A) = 



2V 



(3(2irl) 4 



l_ C + ln(-^-) ) ln(AAr d ). (23) 
Dr d 



in overall agreement with the result of Gefen and At- 
land p3| for two-dimensional non-diffusive systems. 

8. In models of the QCD vacuum such as the instan- 
ton liquid model, a chiral transition may be envisioned 
by varying the density of instantons in fixed V. The 
low density phase is characterized by instanton clusters 
where the quarks and antiquarks are localized over the 
cluster sizes. At high instanton density, the quarks are 
delocalized and chiral symmetry is spontaneously broken. 
These transitions are actually quite generic and may be 
triggered by changing the external parameters in QCD, 
including the shape and size of the Euclidean box V. For 
instance, the formation of quark clusters is strongly fla- 
vor dependent, so it is natural to assume that such a 
transition may take place by just varying the number of 
flavors |24|,^5| . In the following, we will assume that these 
transitions are of the metal-insulator type and proceed to 
analyze their signature on the quark spectrum. 

Let £ be the correlation or localization length in 
the cross-over region. For £ < L, the diffusion be- 
comes anomalous with a scale dependent diffusion con- 
stant [p6|j2^]. The scale dependence maybe estimated 
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using renormalization group arguments ]q,p8|. At the 
critical point with £ ~ L and s»l, the level correlation 
asymptotes R(s) ~ s -i-VM) The critical expo- 

nent v sets the rate of divergence of £ in terms of the crit- 
ical conductance o - *, £/Z « |1 — ct;/ct*| _1/ ||. In pq] this 
result was reached by resumming the multi-2-diffuson di- 
agrams and using renormalization group analysis. Their 
arguments carry to QCD by trading 2-diffusons with 2- 
pions (see above). In particular d = 4, j3 = 2, and 
a* w -F*£ 2 - From the 2 + e expansion flsj,p8|l, v = 1/2 and 
the pion decay constant at the mobility edge i 7 "* vanishes 
with universal exponent v. 

At the critical point the normal diffusive regime with 
constant D disappears. Critical disorder is characterized 
by (j* ~ ai The spectral rigidity in the regime : 
cr* < s < er»(Ai/A) with A; = 1/pl 2 , follows from the 
null spectral sum rule through the ballistic contribution 
(E2). Hence, 



S 2 (iV) 



, 1 - C + ln4 N 

/3(2tt) 4 ~* 



XN 



(24) 



The level compressibility \ is weaker than Poisson (x — 
1) and universal. This result is in agreement with the 
original analysis in |^9| (see also |p0[). 

Finally, the asymptotics of 7?(sjfor s ^> 1 in the crit- 
ical regime, implies that the level spacing distribution 
P(s) [[15[ in QCD changes at the mobility edge. The 
result is a modified Wigner surmise that is intermedi- 
ate between Wigner-Dyson and Poisson ^6j,^7 29 1. This 
transition between chaos and integrability suggests that 
in the same regime the quark wavefunctions in fixed V are 
multi-fractal |31f| with a multi-fractal exponent rj = 8\- 

9. The present results may be checked numerically by 
studying quark spectra using lattice QCD simulations, or 
continuum models such as the the instanton liquid model 
or the Nambu-Jona-Lasinio model. They can be readily 
extended to matter since (|l]) is protected by symmetry, 
and supersymmetric gauge theories with flavor p2[ . 
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